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We study a non-commutative deformation of general relativity based on 
spectral invariants of a partial differential operator acting on sections of a vec- 
tor bundle over a smooth manifold. We compute the first non-commutative 
corrections to Einstein equations in the weak deformation limit and analyze the 
spectrum of the theory. Related topics are discussed as well. 
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1 Introduction 

The basic concept of general relativity is the space-time, which is a smooth manifold 
with a globally hyperbolic metric g^^v As it is well known that the dynamics of the 
metric is described by the Einstein-Hilbert action functional (with cosmological con- 
stant) 

S ^ dxg'l\R-2K), (1.1) 

M 

where dx denotes the standard Lebesgue measure on M, g = Atlg^y, R is the scalar 
curvature of the metric, G is the Newtonian gravitational constant and A is the cosmo- 
logical constant. 

It is a rather old idea that the metric can be described by the properties of light 
propagating in the spacetime (that is, by studying the properties of the wave equation) 
and the action of general relativity can be 'induced' by spectral invariants of a second 
order partial differential operator of Laplace type (for more details, see (9]). In a series 
of recent papers [9 8 10| a new theory of gravity has been developed, called Matrix 
Gravity, in which the main idea is to describe the gravitational field by studying the 
propagation of fields with some internal structure rather than light. This idea follows 
from the fact that at a more fundamental level (short distances or high energies) the role 
of the photon could be played by a multiplet of gauge fields. In the papers ||9] (TO] the 
action was constructed by directly generalizing the language of Riemannian geometry 
to the non-commutative case. In [16] we studied the non-commutative limit of the 
theory proposed in ||9][T0l. It turned out that this deformation was non-unique. That 
is why, in ['81 a new model, called Spectral Matrix Gravity, was developed, which is 
based on spectral invariants of non-Laplace type partial differential operators. In the 
present paper we study the weak deformation limit in Spectral Matrix Gravity. 

We would like to stress, at this point, that our approach to deform General Rel- 
ativity is different from the more standard approach of noncommutative extension of 
gravity on noncommutative spaces. In flat space one usually introduces noncommuta- 
tive coordinates satisfying the commutation relations 

[x^,x']^e>'\ (1.2) 

where 6**^^ is a real constant anti-symmetric matrix, and one replaces the standard alge- 
bra of functions with the non-commutative algebra with the Moyal star product 

fix) * ^(;c) = exp(^0^''A|.j/(^+3,)^(, + ,)^^^^. (1.3) 

An extensive review of different realizations of gravity in the framework of noncom- 
mutative geometry, especially in connection with string theory, can be found in Il23l . 
An extension of the star product and noncommutativity from flat to curved space-time 
can be found in [T,T9'|. 

We list below the most relevant differences between the two approaches, a more 
detailed and extensive discussion can be found in ll9l [T0lfm[T6l . 
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The biggest problem with the curved manifolds is the nature of the object ff'^ . All 
these models are defined, strictly speaking, only in perturbation theory in the deforma- 
tion parameter, that is, one takes (F^ as ?l formal parameter and considers /orma/ power 
series in 0'"'. In our approach to Matrix Gravity the deformation parameters are not 
formal and the theory is defined for all finite values of the deformation parameter. 

In the standard non-commutative approach the coordinates themselves are non- 
commutative. This condition raises the questions of whether the space-time has the 
structure of a manifold, and how one can define analysis on such spaces. Moreover, one 
needs a way to relate the non-commutative coordinates with the usual (commutative) 
coordinates. In Matrix Gravity we do not have non-commutative coordinates. Our 
space-time is a proper smooth manifold with the standard analysis defined on it. 

In non-commutative geometry approach the deformation parameter ff^''' is non- 
dynamical, therefore there are no dynamical equations for it. This poses the question 
of what kind of physical, or mathematical, conditions can be used in order to deter- 
mine it. In addition, in many models of non-commutative gravity (as in |2)) is 
a non-tensorial object which makes it dependent on the choice of the system of co- 
ordinates. This feature leads the theory to be not invariant under the usual group of 
diffeomorphisms. In ||2| the authors construct all the relevant geometric quantities like 
connection, curvature, etc. in terms of the non-commutative deformation parameter. 
In this framework they obtain an expansion of these quantities and of the action up 
to second order. In their approach the object ff^^' is a constant anti-symmetric matrix 
(not a tensor). This violates the usual diffeomorphism invariance, Lorentz invariance, 
etc. for which there exist very strict experimental bounds. The main result in |2J is 
the derivation of the deformed Einstein equations. The zero-order part (Einstein) is 
diffeomorphism-invariant, and the corrections (quadratic in theta) are not. Therefore 
the theory contains some preferred system of coordinates and its whole content depends 
on it. Of course, the theory needs to justify the choice of such system of coordinates. 

In our approach, we do not introduce any non-tensorial objects. As a result our 
theory is diffeomorphism-invariant. So, there are no problems related to the violation 
of Lorenz invariance, etc. and there are no preferred systems of coordinates. Moreover, 
the non-commutative part of the metric in our approach is dynamical. We have non- 
commutative Einstein equations for it. The goal of this paper is, in particular, to derive 
these dynamical equations in the perturbation theory. 

In IIT9 1 the authors assume that ff^^' is a covariantly constant tensor. But then, there 
are strict algebraic constraints on the Riemann curvature tensor of the commutative 
metric (obtained by a commutator of second covariant derivatives). In Matrix Gravity 
such algebraic constraints are absent — the commutative metric is arbitrary. 

In the usual approach of non-commutative geometry (as in, for example, |2J), when 
one defines the affine connection, the covariant derivative, the curvature and the torsion 
the ordering of factors is not unique. There is no natural reason why one should prefer 
one ordering over the other. That is, the connection coefficients can be placed on the 
left, or on the right, (or one could symmetrize over these two possibilities) from the 
object of differentiation. Another aspect of the ordering problem is the fact that there 
is no unique way to raise and lower indices. One can act with the metric from the left 
or from the right. Spectral Matrix Gravity, instead, is pretty much unique. There is 
no need to define the affine connection, the covariant derivative, the curvature and the 
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torsion. There is no ordering problem. 

The definition of a "measure", in standard non-commutative geometry, as a star 
determinant (as in |2 1) does not guarantee its positivity. It only guarantees the positivity 
in the zero order of the perturbation theory. With our definition, the measure is positive 
even in strongly non-commutative regime. 

Moreover, the Moyal star product is non-local which makes the whole theory non- 
local with possible unitarity problems. In Spectral Matrix Gravity approach the action 
functional is a usual local functional of sigma-model type (like General Relativity, but 
with additional non-commutative degrees of freedom). There may be problems with 
the renormalizability (which requires further study) but not with unitarity. 

We consider the space-time as a smooth manifold in which the gravitational field is 
described by propagation of gauge fields rather than light. Therefore one is naturally 
led to consider non-Laplace type partial differential operator for which the leading sym- 
bol is now a matrix. One can evaluate the heat kernel asymptotic coefficients for such 
operator and construct, from the first two, the action for Matrix Gravity as in (II. 10b . 
Similar calculations have been performed in noncommutative geometry regarding heat 
kernel asymptotics expansion. In 11211 l22l the author evaluates the relevant geometric 
quantities from an approximate power expansion of the trace of the heat kernel for a 
Laplace operator on a compact fuzzy space. In [24| the author studies the quantiza- 
tion of noncommutative gravity in two dimensions by considering a noncommutative 
deformation (using the Moyal product) of the Jackiw-Teitelboim model for gravity. In 
this case the path integral can be evaluated exactly and the operator for the quantum 
fluctuations can be found. Once the operator is known one can study the first two heat 
kernel asymptotic coefficients and obtain information about the conformal anomaly 
and the Polyakov action. Our analysis is rather different because we consider a stan- 
dard smooth manifold manifold and a non-Laplace operator with matrix valued leading 
symbol without considering Moyal product. 

For the sake of completeness we describe, briefly, the framework of Matrix Gravity 
(a more detailed discussion can be found in |l9] [TO] |8]). Let us consider a smooth 
compact orientable n-dimensional manifold M without boundary. Let y[w] be a A^- 
dimensional vector bundle over M of densities of weight cj and C°°(y[w]) be the space 
of smooth sections of the vector bundle V{<x>\. Locally, a section of the vector bundle 
V is described by the A^- vector i/A. Let End(y) be the bundle of endomorphisms of 
the vector bundle V; a section Q of the endomorphism bundle is represented locally 
by the matrix Q^b- Next, we assume that the vector bundle V{<x>\ is endowed with an 
Hermitian metric G of weight (1 - 2aj) defining the fiber inner product - <p^Gi//. 

This enables one to define the L^-product of two sections of V by 



The completion of the space C°°(y[<jj]) in the corresponding norm defines the Hilbert 
space L^{V[(L)]). 

Now, let us consider a second order partial differential operator 




(1.4) 



M 



L : C°°(V[lo]) C°°(yM) , 
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acting on smooth sections of the vector bundle V[<x>\, with the endomorphism- valued 
coefficient functions 

L = -ai'^xjd^.dy + b^{x)df, + c(x) . (1.5) 
The leading symbol of the operator L is defined as 

crdL;x,^)^a^'\x)^^^y, (1.6) 

where is a cotangent vector at the point x. 

We consider two cases: the Euclidean case when the operator L is elliptic, and 
the pseudo-Euclidean case when the operator L is hyperbolic. These cases are related 
by an analytical continuation and, for the sake of simplicity, we will restrict in our 
calculation to the elliptic case. The formulas for the hyperbolic case look exactly the 
same (for more details, see 19][l0l[8l). The operator L is elliptic if the leading symbol 
o-l{L; X, ^) is non-degenerate for any x and any ^ + 0. The operator L is of Laplace 
type if the leading symbol is essentially scalar, that is, 

cr(L-x,0-lg''\x)^,^, , (1.7) 

where I is the identity matrix. 

The geometric invariants we need can be derived from the spectrum of the operator 
L. The L^-trace of the heat semigroup, exp(-fL), usually called the heat trace, contains 
all the information about the spectrum of the operator L. It is well known that the 
asymptotic expansion of the heat trace as f — > 0^ has the form fTTl 

oo 

Tr^. exp(-fL) ~ {An)-'i ^ f*-Ut , (1.8) 

where the coefficients Ak are called the global heat kernel coefficients, which are, of 
course, spectral invariants of the operator L. They have the form 

Ak-^dx tryak , (1.9) 
M 

where try is the fiber trace on the vector bundle V and the coefficients a^ are some 
endomorphism-valued densities of weight one, called the local heat kernel coefficients. 
It this the coefficients Aq and Ai that we study in this paper 

The action of Spectral Matrix Gravity proposed in IfTOl is a linear combination of 
the heat kernel coefficients Aq and Aj, more precisely, 

^ [6Ai-2AAo] . (1.10) 



16;rGA^ 



We would like to point out here that the above action can be also thought of as a 
particular case of the Spectral Action Principle introduced in the framework of non- 
commutative geometry in 1 1 3] and fl^ . 

For the Laplace operator, L - -A, the heat kernel coefficients are well known 

1 

aQ=I, ai = -RI. (1.11) 
6 
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In this case, the action of Spectral Matrix Gravity reduces to the standard Einstein- 
Hilbert action ( 11.11 ). 

We would like to stress, here, that we are interested, in this paper, in a much more 
complicated general case of an arbitrary non-Laplace type operator (with a non-scalar 
leading symbol). In this case there is no preferred Riemannian metric and the whole 
language of Riemannian geometry is not very helpful in computing the heat kernel 
asymptotics. That is why, until now, there are no explicit general formulas for the 
coefficient Ai. A class of so-called natural non-Laplace type operators was studied in 
EE] where this coefficient was computed explicitly. 

The main goal of this paper is to study the action of Spectral Matrix Gravity in 
ffie weak deformation limit and to describe the corresponding corrections to Einstein 
equations. For the sake of completeness we will briefly describe, in the next section, 
the properties of some geometric two-points quantities that are widely used in the cal- 
culation of the heat kernel asymptotics. 

2 Geometric Background 

We introduce now some geometric two-point quantities following mainly |T, "Sl (see 
also, ||20l [T5l |5] [3] lU). Let x' be a fixed point on the manifold M. There always 
exists a sufficiently small neighborhood of the point x' such that every point x in this 
neighborhood can be connected to the point x' by a unique geodesic. In the following, 
then, we will restrict ourselves to this neighborhood. Moreover we will denote tensor 
indices at point x' by primed letters. The non-primed and primed indices will be raised 
and lowered by the metric at the points x and x' respectively. Also, the derivative with 
respect to the coordinates of the point x' will be denoted by primed indices as well. 
We will use the standard notation of square brackets to denote the coincidence limit of 
two-point functions, more precisely, for any function of x and x' we define 

[/](x)= lim/(x,x'). (2.1) 

First of all, we define the world function cr{x,x') as one half of the square of the 
length of the geodesic between the points x and x' . Next, we define the first and second 
derivatives of the world function 

o-^ = V^cr, cr^, = V^'Cr , (2.2) 

= VvVo- if'y = VvVV , (2.3) 
and the Van Vleck-Morette determinant 

A(x, x') = e^^^^'^"'> = g-Hx) det[-V^Vv'0-(x, x')]g"^ (x') , (2.4) 

where ( - log A'^^. Now, let be a section of a vector bundle V, be a connection 
on the vector bundle V and "R^y be the curvature of this connection. Next, let fix, x') 
be the parallel displacement operator of the section (p along the geodesic from x' to x; 
obviously. 
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Finally, we define the first derivative of the operator of parallel transport 

=!p-iV,!P. (2.6) 
One can show that these functions satisfy the equations |I3]|5| 

(T = \(T,(T^ = \(^'(r,, , (2.7) 

f/'^cr" = C7^ , Tf^cr'^d^', Tf'ycr^, = cr, , (2.8) 

cT^V^f = i(«-^V- (2-9) 

cr^V^!P = <7^'£^ = , (2.10) 

and the boundary conditions 

M = K] = K'] = 0, (2.11) 

[fv]=5^, [rf'y]^-S^y. (2.12) 

[A] = l, [^] = 0, (2.13) 

[!P]=I. (2.14) 

The coincidence limits of higher derivatives of the two-point functions introduced 
above are expressed in terms of the curvature. In particular, the ones that we will need 
below are EElllSlIll 



[C,^,] = (2.15) 

K^v] = ^Rt^v, (2.16) 

[rf'v;a] = 0, (2.17) 

[rf'y.a/s] = -^(R''val}+R''ayis) , (2.18) 

[Ey] = (2.19) 
1 
2 



[Ey-^] = --K^, (2.20) 



where R^yp^ is the Riemann tensor and R^jy = R^jjay is the Ricci tensor. 



3 Heat Kernel 

We will use a method for the calculation of the heat kernel developed in UISI, which 
is based on a covariant Fourier transform proposed in |l3]|5l. The heat kernel for the 
operator L is the kernel of the heat semigroup, that is. 



U(t\x, x') = exp{-tL)r(x, x')6(x, x) , 



(3.1) 
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where 5{x, x') is the deha-function (in the density form). Now, following ||3] we define 
the CO variant Fourier transform as follows. First, we define 

(5(x,jc') = A5(x,x') f exp{(^/,-cr^'(x,jc')) . (3.2) 



Then, by using this representation and the equation 

exp{-fL) / = / exp{-f(/-iL/)) , 

we obtain 

U{t\x,x')^ ^^x,x')'P{x,x') f -^^&x'p{i^i,,(j^'{x,x'))<b{t\k,x,x) , (3.3) 

where 

(^{t\k, X, x') = exp(-fA) ■ I , (3.4) 

A = e-''^^'"''' ^-^-L^^'Pe•^^^"^' . (3.5) 

By using the coincidence limits of the two-point functions we obtain the heat kernel 
diagonal 

U(t\x,x)^ r mKx,x). (3.6) 
J (2n)" 

Now, we represent the operator L in the form 

L = -p-'V^pa'^'pVyp-' + Q , (3.7) 

where a^^ is a matrix-valued symmetric tensor of type (2, 0), p is a matrix-valued den- 
sity of weight 1/2 and Q is a matrix-valued function. By substituting the operator (|3.7| l 
in equation ( 13. 5t , we get 

A = -e-'^''''^' A-ip-^'p-^V^'Ppa''''pp-^VyPp-^ A^-e'^^-"^' + Q , (3.8) 

where a'"' = ' ai'^'P, p - 'P'^p'P md Q - P'^QP. We rewrite this operator in a more 
convenient form 

A = -l^a^X + Q , (3.9) 

where 

= e-'f'''""'''A-^p-i!p-iv^!PpA^e'f""^' , (3.10) 
It is useful to introduce, now, two quantities 

Cv = -P;vp' and Cv--p'^P;v- (3.11) 

Then we get 

Xy ^Vy+Cy+^;y+Ey + l^^'lf'y , (3.12) 



G. Fucci and I. G. Avramidi: Spectral Matrix Gravity 



8 



X^^N^-C^+(.,^+E^ + i^,,r]\, (3.13) 

and 

A = -(V^ - + + + l^p,rf'^,ra^'^V, + C, + C,v + E, + i^p^rf' y) + Q . (3.14) 
Finally, a straightforward calculation gives 

A^H + K + L. (3.15) 



Here 



where 



H = ,if\a^'\ (3.16) 

K = -!^p,(B'''''V, + G^') , (3.17) 
£. = -Dpa^'^'Dy + Q, (3.18) 



+ Epa^'if'y + i/'yO^'Ef, + l^.^a^^if'y , (3.19) 

©V = Vy+Jly^Vy + Cy + (;y + Ey, (3.20) 



with Cv, C^, { and fi^ defined in dTTTT i, ( |24] | and dZg] ). 
More exphcitly we can also write that 

X = -a/^-V^V, + J/f + Z , (3.21) 

where 

J/" = -a""; , + Cva/^" - a'^'Cv - 2a'"'^; , - a'^''£v - £va''\ (3.22) 

- EpCf'Ey - ^; .a^^E^ - ^; yE ^ci''' + ^ . (3.23) 

Thus, by using the eq. ( 13.151 1 we obtain 

U{t\x,x)^ r JL.e-''^-^-^^-\ , (3.24) 
J (In)" \x=x' 

R" 

which, by scaling the integration variable ^ — > r^^, takes the form 

t/(f|jc, x) = (47rf)"^ r ^ exp(-// - ^ftK - t£) ■ l| . (3.25) 
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It is convenient, to rewrite this equation as 

U(t\x,x) = i4-7:tr-^ f ^ e-^^^'expi-H- yrtK-t£)-l\ , (3.26) 

R" 

Where l^p = gi^'^^^, and 

H = H- . (3.27) 

In order to evaluate the first three coefficients of the asymptotic expansion of ( |3.26l l 
as f ^ we use the Volterra series for the exponent of a sum of two non-commuting 
operators 



We obtain 

exp(-H - yftK - t£) ^ e-" - Vf Q + f^F + 0(t^ , (3.29) 



where 



1 





and 

1 T2 1 



^ = J dTie-^^-'''''" Ke-^'" , (3.30) 



T = . (3.31) 



We are only interested in the terms oq and ai of the heat kernel expansion, namely 

the terms of zero order and linear in the parameter f. These terms can be written, 
respectively, as 

flo - g-ao, (3.32) 

fli = , (3.33) 



where 



ao 



a\ — 



f^g'-^ e^lf'' exp{-//) ■ l| , (3.34) 
f ^ g--2 e^^^'*^' ■ l| . (3.35) 

J 712 l.t=.v' 



The term of the heat kernel expansion vanishes identically. This happens because 
the heat kernel coefficients are defined as ^-integrals over the whole R" and the term Q 
in ( 13.30b is an odd function of ^. 
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4 Evaluation of the Heat Kernel Coefficients 
4.1 Local Coefficient uq 

We will evaluate the heat kernel coefficients Aq and A\ using the perturbation theory. 
The main idea is to introduce a small deformation parameter A and evaluate the non- 
commutative corrections to the action of Spectral matrix Gravity. For this purpose we 
write the matrix a^y as 

a"' ^g^n + Ah''^ (4.1) 

where /z'"' is a traceless matrix-valued tensor field (a non-commutative perturbation of 
the Riemannian metric), satisfying 

txyh^'^ = . (4.2) 
Furthermore, we parametrize the matrix-valued density p introduced in ( 13.71 ) as 

p = gie'''e^'' . (4.3) 

Here cr is a traceless matrix-valued scalar field and is a scalar field. {Do not con- 
fuse it with the world function introduced in the previous sections!) Finally, we also 
decompose the endomorphism Q, 

Q^ql + A®, (4.4) 

where © is a traceless matrix-valued scalar field. 

Now we expand all the quantities in powers of A. On doing so the matrix p and its 
inverse read 



p = g-* e1'[\+ Aa-+^a^'^ + 0{A^) ^ 



p-^ = g-i e-'^{l-Ao-+ ^(T^ ) + O(A^) , (4.5) 



and its derivative is 



A^ 

^o-; V + -^(o"; vO" -H crcr- v) 



+ e'^<P,y\l+Ao-+YO-\ + 0{A'). (4.6) 



From the last two expressions one can easily evaluate the operators Cy and Cy obtaining 
explicitly 

Cy = -P, yp-' ^-Cf>,y- AH, V + ^ [^; V, + 0{A^) , (4.7) 



Cy = -p-'p; y^-Cf>,y-Aa;y-— [^; y, (X] + O(A^) , (4.8) 

where cr - V^^crV. 
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The operators H, K and £. introduced above in (13.161 1, (13.171 1 and ( 13.181 1 depend on 
the deformation parameter A as well. By expanding them in terms of the deformation 
parameter we get 

H = Hq + AHi , 

K = /To + ^Ki + A^K2 + O(A^) , 

£ = U + A£,+A^£2 + 0{A^), (4.9) 
where (after defining = 'P-^h^"''P and = !P '0!P) 

Ho = ^a'^Hrj-',if\g'"'-8"'^'), (4.10) 

H, = (4.11) 

Kq = -i^„-(2;7''VV, + 2/^g'^''^;, + / ^^'^ + £f/^ + 7,«>'^), (4.12) 

+ J]"',E,hn, (4.13) 

^2 = -'^.'['7";([^^;v,/in + [c^'^ (?■])] , (4.14) 

Xo = -V2-(2^'''+2£OV^+0.^'''+0;^f "-^.^•^ 

- ^■''f;^-£/'^-£''£^-2^;^£''+^, (4.15) 

Xi = Vv-[/i''^v + 2^; v/i'"' + (F"£v + £vO]V^ + ^;/'^ 

+ F''; ^0; y + /^V; + 20' "OT; ^ + 0; ^h""^, y " /i''"; y 

- F^^; - ^; ^/l'^''^; , - ,Ey - h'''' E y, , + [E„&^''] 

+ (f>-^,[Ey,h^"']-Ef,h^"'Ey-E^h^''Cy-(,,h^"'Ey + @, (4.16) 



X2 = W\ y] -[<T''',d-]W, a] +h^r^^a,y + h^^ 



+ E^Ji^^'a^y-a-.^^'Ey . (4.17) 

In the framework of perturbation theory we write, then, the coefficients ciq and 5i 
of the heat kernel expansion in (I3.34l i and ( 13.351 1 in terms of the deformation parameter 
A, namely 

M^) = af + Aa^^^ + A^af + O(A^) 

axiA) = af^ + Aa-^^ + A^af + 0{A^) . (4.18) 

By using the explicit formulas obtained in (14.10b through (14.17b . we will be able to 
evaluate all the coefficients of the Taylor expansions in (14.18b . 
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Next, we introduce a notation that will be useful in the following calculations. Let 
/ be a function of ^. We define the Gaussian average of the function / as 



</> 



j n'i • 



(4.19) 



The Gaussian averages of the polynomials is well known 

(2n)! 



(4.20) 



where the parentheses ( ) denote the symmetrization over all the included indices. 

For the coefficient of order zero of the heat kernel expansion we consider the first 
equation in ( 14.181 ). From the formula (13.341 ). it is easy to see that the only non-vanishing 
contribution to oo is 



ao ^ {\\ - AH, + —H{ 



(4.21) 



By using the equations ( 14.10b . ( 14.1 11 1 and after taking the coincidence limit we obtain 
the expression 



I - Ah""^,^, + ^A^f^'hP^^.^y^p^^ I ) + OiA') 



and then, by performing the Gaussian averages, we get 

A A^ 

5o = 1 - + —ih^ + Ih'^'h^y) + O(A^) 
2 o 



(4.22) 



(4.23) 



where h = gfivh^^- In order to evaluate the global coefficient Aq, given in il.9i . we need 
the trace of ( 14.231 1. Since h'''^ is traceless we immediately obtain 



(4.24) 



tryflo = try 1 1 + —h^ + — I + OiA') . 



4.2 Local Coefficient ai 



Now we evaluate the coefficient a,. By using the expressions ( |3.35t . ( 13.311 ) and ( 13.181 1 
we have 



tryfli 



1 T2 

(try [-e-"Q]) + Itrv J dr2 j dne-^'-'^-^^Ke-^'^--^^^" Ke'^'" | 



try 



1 

/ 



dT,e-'^-'"^"£).,a''''Dye-''" 



(4.25) 



G. Fucci and I. G. Avramidi: Spectral Matrix Gravity 



13 



where the first term of the expression (14.251 1 has been obtained by simply using the 
cyclic property of the trace. 

In the following we will evaluate the terms in ( 14.251 ) separately. We start with the 
simplest of them, namely the one involving the endomorphism Q. By using the Taylor 
expansion m Aof H in ( 14.91 ) and the coincidence limits ( 12.121 ) we obtain 



(4.26) 



(try {-e-"Q\)\^^^^ = {txy {-Q + A^.^^^'Q - L^^^^^^^^h^^hP^ 

By expanding Q as in (|4.4| l and by performing the Gaussian averages we obtain 

(try {-e-"Q\) I = -Nq + Ahvy i^-h® - ^h^q - ^h'^'h^yq'^ + 0{A') , (4.27) 



where we used the property ( 14.21) . 

For the second term in equation (14.251) we get, by using the definition (13.171 1. 



try 



1 T2 



= - try 



dTie-^^-^''>"[BP'''Vye-^^'-^'^"B'^'''Vf,e-''" 



^p'^o-' J dT2 J 


+BP'''Vye-^^'--''^"G'^'e-''" + G^' e-^'''-^'>" B'^'i'V^e-^'" 



+GP'e-^^'-^'^"G'''e-^'") ^ 



(4.28) 



It is straightforward to notice that in the last expression we need to compute first and 
second derivatives of the exponentials containing the operator H. These derivatives are 
computed by using integral representations, i.e. for the first derivative we have ifTOll 



tH 



where 



r 

J ds e-'"H,^ 



(4.29) 



(4.30) 



This last integral can be evaluated by referring to the following formula and by inte- 
grating over s IfTOll 



k=0 



k\ 



By expanding (14.301 1 in s, up to the second order in A, we obtain 

1 



(4.31) 



(4.32) 
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This last expression can be obtained by recalling that the coincidence limit for H and 
its derivatives is of order A without the zeroth order term (see Appendix A). 
For the second derivative we write [,10J 



T 



T S2 

+ JdsiJ dsi[e-'-''--"^"H,ye-'""H,^e-^' 



-(t-S2)H 



We can express this formula in the same form as (14.32b . i.e 

1 ~ ~, 1 



(4.33) 



(4.34) 



Now that we have the expressions ( 14.291 ) through (14.34b . we can substitute them 
in (14.28b and we can expand the remaining exponentials in t up to orders A^. After 
the expansion of the exponentials and after taking the coincidence limit, we have to 
evaluate the double integrals of polynomials in ti and T2 which will yield the numerical 
coefficients for the various terms in ( 14.28b . The most general double integral that we 
need to evaluate is the following 

1 T2 

h{a,P, y,d)^ J dT2 J dnT'lil - Tif{T2 - T{f(\ -T2+ TlY , (4.35) 


with (a,/5, 7, 6) being positive integers so that the integral is well defined. The solution 
in closed form of ( 14.35b is (see Appendix B) 



r(i + Q')r(i + yS)r(i + r)r(2 + a+i3 + d) 

r(3 + a + yS + y + 6)r(2 + a + /3) 



(4.36) 



where r(jic) is the Euler gamma function. 

By using the technical details described above and the coincidence limits for the 
various terms in ( 14. 28b (see Appendix A), we obtain 



try 



j ^T2 J dTie-^'-'^^"Ke-^''---'^"Ke-''" 



= tiv i^o) + ^ir-v (Qi) + ^^try (Qj) + O(A^) 



(4.37) 
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where 



1 

-R 

6 

, 

1 



(4.38) 
(4.39) 



1 



1 



Q, = h-Jri" + -h^'". „h. „ - -/z'^". nhuv' P h\i^.„ + —hti"^ 



16 



12 



12 



We can finally evaluate the last term in equation (|4.25l l. By using the definitions 
(13.181) and (13.201) we can write that 



try 



/- I) 



1 

-- (try j t/Tie-<'-^'^^(fl^^^V, + fl'^'V^V, + a'''',^Jly 



(4.41) 



In order to evaluate this term we use the derivatives in ( |4.32| l and ( |4.34t and we 
expand the remaining exponentials of ^ in t up to terms in A^. During the calculation 
the numerical coefficients of the various terms can be evaluated by referring to the 
following general integral 



/2(a,;e) = j dTlT^^d - T^f , 



(4.42) 



where {a,/3) are positive integers, and for which the general solution in closed form is 
(see Appendix B) 

r(i + Q')r(i +5) 

The explicit form of ( 14.41b can be obtained with the help of (14.431 ) and the coinci- 
dence limits in Appendix A. After a straightforward calculation one gets 







= try (Ho) + ^trv (Si) + £-t[v (H2) + 0{A^) , 



(4.44) 
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where 



-d,. .'f - d,. .rf,-f (4.45) 



-0 - -(P;tJ -9;tiV 



^1 



= , (4.46) 
4 z o 4 



+ hO-;^(f>'^' + -hh^"'(f,, ^0; V - Cr. ^Cr- - V; ;,<^; y . (4.47) 



In the notation of equation ( 14.18b we can write, now, the different contributions, in 
increasing order of A, to the coefficient ai. In more details, by using the results ( I4.27l i. 
(14.381 ). ( 14.451 ) and recalling that 

trvai(^) = tryaf^ + itrya^/^ + yof + O(A^) , 

we get 

af^ ^^R-<p,^-'^-cl>,^<p-'^-q. (4.48) 

Moreover, by using ( |4.39l l and ( 14.461 ) we obtain 

fl^i" = 0. (4.49) 

Finally, by combining the results in (14.271 ). (14.40b and ( 14.471 ). we have the following 
expression for the term of order inai, i.e. 



a 



(2) 



^h>'lZf,y<f,,p4>'P + j^h;^h''' - ^hhP'Rf^y - ^-h%h'PR^,y + ^/Z^/? 



+ —hf^huyR + -h& . (4.50) 
24 2 



5 Construction of the Action 



In order to write the action of Spectral matrix Gravity, we need to evaluate the global 
heat kernel coefficients Aq and A; . As we already mentioned above, the coefficients Aj. 
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are expressed in terms of integrals of the local heat kernel coefficients aa (which are 
densities) or the coefficients a^ (which are scalars). Namely, 



M 

By using the equation (14.241 1. we get 

i2 



dx try 



(5.1) 



(5.2) 



Now we use the equations ( 14.181 1 and ( |4.48t - (|4.50t to compute the coefficient Ai. 
By integrating by parts and by noticing that the trace of a commutator of any two 
matrices vanishes, up to terms of order A^, we obtain 



J dx g'^ tivl 

M 



q- 



1 7 



o-'''cr.„ + \h& 



12 6 12 ''^ ' 12 48 



—h^^'h R - —h h'^ - —h^"'''Ph 



0{A') 



(5.3) 



The invariant action functional is written as Unear combination of the coefficients 
Aq and Ai as shown in (ll.lOl l 



M 

— try^ - 6o-'''cr.^ - 1\vP(T,y, + 3/z0 + 3M''>;^0;v + 6/jCr' 

ih^^h^cj), y - Uo-, yh""^,^ - ^hh^'R^y + ^h'"'hP''R^^py - ^h^hP^R^y 
\h'"';Av + \hPr^yh^p-P + ^h^R + ^h^lz^yR - ^h.^lr" - ^h^^^'^h^y^p 
^-h'q-^-h%^'^' - ^-h\,4>'' - \h''Kyq-\h^^h,y<^,p-P 



A, 



A 



huyh"' 



+ 0{A'). 



(5.4) 
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It is straightforward to show that the action functional that we obtained is invariant 
under the diffeomorphisms and the gauge transformation /i'^" — > UM'^U'^. 

The next task is to find the equations of motion for the fields <t, /j^v and <p by 
varying the action functional. In this way we will explicitly find the noncormnutative 
corrections to Einstein's equations. 

6 The Equations of Motion 

By performing the variation with respect to the field cr, we obtain the equation 

Ana + nh- Ihn^ - Ik' >; v + AW"- y<p. ^ + Ah^^tp- + 0{A^) = . (6.1) 

Here □ is the Laplacian in the Euclidean case and the D'Alambertian in the pseudo- 
EucUdean case. For the matrix- valued field we obtain the equation 

o 



2 ' '63 

(6.2) 

The variation of the action with respect to the scalar field ^ yields 

i2 , 

- 2h " - IhUa + F"; v/j; ^ + /i^"/?; ^.V + 40-; ^.yW + 40-; vF"; ^ 

- hiuh - ^h,^,h"' + + ^h^a4> - hf^ah^y - h^y^plf^-P 

+ 2h^yhi'^' P(p, p + hf^hf.yOcpj + O(A^) . (6.3) 

The equation of motion for the field g''^ can be written in the following form 

1 /l^ 
R^" - ^g^^'R + J^g"" = T""" + —tivx/^"' . (6.4) 

Here the tensor J'^" is 

T"" = 60- ^cf)' - Sg^^cf), p<p- P - Sqg^" , (6.5) 

which represents the stress-energy tensor for a massless scalar field. The tensor JH'"' 
represents, instead, the stress-energy tensor for the fields h^^ and o". 
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The equation (I6.4l l is the main result of our paper As we can see, the new fields 
of our model, /i'''' and cr, contribute to modify the standard Einstein equations. More 
precisely they contribute to an additional term in the stress-energy tensor 

The tensor jz/'"' can be written as the sum of six terms: 

=5/^^ = + ^i^^ + + + + . (6.6) 

In the first term we have only derivatives of the field cr 

j^^^^ = 6cr' <^o-' - Sg^V; po-' P - 6ho-' <^0; + 3o-' <^/i' ''^ + 6// V' ^<^; p 

+ 3cr;p''^/i'"' + 3/" J/!o-''^f p - ^o-'Ph,p - lo-.p^.^hP"^ . (6.7) 
The second term only contains derivatives of the scalar field 0, namely 

= 3//>; p0; , + i^^^/l j + 3^; , j/l'^^ p//" " ^g^"' h, 

- ^-{cp, pcf>' P + cf>,p- P) ^Ih^Ji^^'- + hh>'^ + i/z V + ^g'^h^f^h,,^"^ 
+ 3hP''h'"'4>,p^ + ^(0- "^^^ + ^- ">) Jv'^'^- + 1/^2 j . (6.8) 
The third term only contains second derivatives of the matrix-valued tensor field h^^, 

= h^'^, (.^j/l")- + i/;'-; .^/Z" - ^-hh"^'^- , + ^-g'^'hhP'^, p^ 

+ i/z'^/z")^; p,. - ^hP'^h'^r^ p^ - /,p,: ClPl/z'^)- + /xf-p/^f 

(6.9) 

The fourth term contains only find first derivatives of /i'^'', namely 

<; = —hP'^^ph'"; . + \lv''h>^',p - \hP'^, (p/z'''^ + \g>^^h, ^hP'^,p 

+ hi'^^,p0p, ^ + hiP'^' ^hp^' - ha-p' '"h^'^""' p + hp^"' i^i/ip''). r 

z 4 z o 

(6.10) 

The fifth coeflicient contains only first and second derivatives of h 

^ ^ 4 4 o 

- \g''{\^hh,pP-h,pJiP''+^-h'Ph,y (6.11) 
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The last term, -^^gj, does not contain any derivative of N^^, namely 



Ifivhl'V + hl^^h^P + gl^v i^hh<^fi + h^h^P) + h^g^^g'^f' 



1 



+ hP'^hp^g'^^^g 



+ -h'"'hP^R^ 



-g'"'hP''h'"^R, 



4 

+ ^(^ - 6?) [g'^hp^hP'^ + Ahp^Ph)'^P) + ^(R - 6q) [g^^'h^ + Ahh^''') 



(6.12) 



The dynamics described by the equations (16.1b . ( 16.21 ). ( I6.3l l and ( I6.4l i can be studied 
by using an iterative method. Let us write the solution for the background fields (p and 
g^^'^ as Taylor expansion in the deformation parameter A as follows 



= g'J + Ag'^' + AY^ + OiA'). 



(6.13) 



By substituting these expressions in equations (16.31 ) and ( 16.41 ) we obtain, for the terms 
of order /l", the dynamical equations 



□0 



. 



(6.14) 



As we can see from the last equations the term g^'' is nothing but the solution of the 
ordinary Einstein equation in vacuum with cosmological constant. By substituting the 
solutions to ( 16.141 ) back into the equations of motion for the fields cr and we get 
equations of the form 



<b2(g^\cf,oW'' 



O(A^) 



(6.15) 



where <l>i(g^'', 0o) and <l>2(g^'^, 0o) are linear second order partial differential operators. 

By iterating this process we can, in principle, find the solution to our dynamical 
equations in form of a Taylor series in A. 



7 Spectrum of Matrix Gravity on De Sitter Space 

The action for matrix gravity obtained in the previous section is a functional of the 
fields (p, cr, h^"^ and g^'^. The dynamics is described by a system of non-linear partial 
differential equations coupled with each other. We analyze, now, the dynamics of the 
theory. For simplicity we will set, from now on, Q - 0. This particular value for the 
matrix-valued scalar Q will not affect our analysis. 
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As already mentioned above, from the equation of motion ( 16.31 ) for the field (fi, we 
can see that to the zeroth order in the deformation parameter A the field satisfies the 
following equation 

□0 = 0. (7.1) 

As it is well known, the solution of the last equation represents a wave propagating 
in the whole space. Since we require that (fi vanishes at infinity, the only solution is 
4> - 0{A?) in the whole space. 

With this solution for the field 0, the matrix-valued function p defined in ( |4.3l l 
becomes 

p^e^'^ . (7.2) 

A deeper analysis shows that the matrix-valued scalar field cr is not an independent 
field. Following ||9][l0l the general form of p can be written as 

p = , (7.3) 

where 



= -^e/.,...^„ ev,...v„ fl''"'' ■ ■ ■ fl^'"""' . (7.4) 

By using the decomposition (14. Il l of /z'^*' in equation ( 17.4b we get the following formula, 
up to the term linear in the deformation parameter A, 



a-^-^h+'-^h^yh'''' + 0{A^) . (7.5) 



We can write down, now, the action by imposing the constraints (p - 0{A^) and ( I7.5l l. 
The final result is the following 



A 



2 



16.G I — I -.^''^;.-^-2A+^tr. 



- g^"'h h 



M 



)} + O(A^) . (7.6) 



The action depends, now, only on the independent tensor fields , W" and the 
scalar field 0. Therefore, we will have only two equations that describe the dynamics 
of the theory. These dynamical equations can be easily derived from the ones given in 
the previous section by imposing the conditions (I7.5l l and = 0{?l-'). 

The action ( 17.61 ) and the equations of motions for the fields evaluated in the previous 
section, assume a simple form on maximally symmetric background geometries. As we 
mentioned in the previous section, the /l" term of the background field is solution 
of the Einstein equations in vacuum with cosmological constant ( I6.14I I. In this section 
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we consider the De Sitter solution to the equation (16.141 ). In this maximally symmetric 
case the Ricci and Riemann tensors take the following form 



R'^vap = , ^ -.A S^agvp - S''/3gya)R and Rf,y = -g^ivR 



n(n — 1) ■ ' ' n 

The De Sitter metric gives a solution of the classical equations provided 

2n 



R 



n-2 



-A 



(7.7) 



(7.8) 



By substituting the expressions in equation (I7.7l i in the action ( 17.61 1. we find a form 
of the action functional valid in De Sitter geometry, namely 

S = j dxg--l^-6<f>'''cl,,^+R-2A 

M 



where the terms /ii and /12 are defined as follows 

«^ - 5n + 1 



/J2 = - 



2n(«- 1) 
n-3 



■R-A, 



n-\ 



-R + 2A. 



(7.9) 

(7.10) 
(7.11) 



It is interesting, at this point of the discussion, to derive explicitly the spectrum 
of the theory. In order to achieve this result we need to decompose the field hf^, in 
its irreducible modes: traceless transverse tensor mode, transverse vector mode, scalar 
mode and trace part. In other words we can write /i^v as 



where the scalar field (f is defined as follows 

tp — h — atff , 
and the fields h^^ and satisfy the conditions 



(7.12) 



^"hj;, = , 



g^^hl;, = , = . 



(7.13) 



We can now substitute the expression ( 17.10b in the action (17.91 1. and evaluate the terms 
separately. Explicitly we obtain 

j dxg-^{^hii'r^hpf,-P^ = jdxg^- - ^ipaip- ^^i^a + ^^atl/ 

M M 

1 , / ,„ 1 / 1 
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M M 

1 / R\ 7 ' 



(7.15) 



/ 

M 



n 



M 

2 , / R\l n+l \ 

-ip(-D+iJ2)Dt(r + 2{^[a + -]\a-ji2 + -R\C 

n ^ \ n I \ n(n - 1) / 



R I 2R 



3R 

□' + 1 ^i2\a + -[ 

n I n \ n— I 



■^2 



Ulf/ 



(7.16) 



and finally 



^ dx g2 h{-u + iii)h = 

M 

^ dx g^-[ip{-u+ nx)ip + 2ip(-n+ iii)n(l/ + il/{-n+ iii)n^il/\ . (7.17) 



By using the decompositions (17.14b through (17.171) we rewrite the action ( 17.91 ) in 
terms of the irreducible modes of /i'"', namely 



/dx ?2 



•'^f ^+7?-2A+ — try 



4 



-^,{-a+fi2)h^'''' 



(»-l)(«-2) 



n _ n(n + 2) 
- r— 



2(n-l) («-l)(n-2) 



3 / 2R\ln-2 \ 



n + 2 In — 2 



An 



-R-2K\nilj 



(7.18) 



It is straightforward to show now that on the mass shell, ( 17.81 1. the terms containing 
the fields and if/ vanish identically. More precisely we obtain the following form for 
the action functional 



1 r , f 4A 

M 



4 («-i)(«-2); 



(« - !)(« - 2) 
4n2 



2«A 



-□- 



(«-l)(«-2) 



+ (9(^''). (7.19) 
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Thus on the mass shell the only remaining fields are the (traceless) matrix-valued trace- 
less transverse tensor h^^ and the (traceless) matrix-valued scalar field (f. This action 
looks exactly the same as in general relativity, the only diff'erence being that the fields 
are matrix-valued and traceless. Therefore, it describes (A^ - 1) spin-2 particles and 
(A^ - 1) spin-0 particles. Note also, that exactly as in general relativity the scalar con- 
formal mode is unstable. 

8 Concluding Remarks 

In this paper we studied a noncommutative deformation of general relativity (called 
Spectral Matrix Gravity) proposed in ifTOl where the noncommutative limit has been 
explicitly evaluated. The approach of the paper fTO^ to construct the action for Matrix 
Gravity differs from the one proposed in |9 8 16|. In the latter the action of our 
model was a straightforward generalization of the Einstein-Hilbert action in which the 
measure and the scalar curvature were matrix-valued quantities. This last approach 
seems to have some intrinsic arbitrariness due to the freedom of choosing the particular 
form of the matrix-valued measure (for a discussion see fT6l ). In order to avoid these 
issues, in flOl the action of matrix gravity was defined as a linear combination of the 
first two global heat kernel coefficients (ll.lOl i of a non-Laplace type partial differential 
operator 

By using the covariant Fourier transform method we were able to evaluate the coef- 
ficients A() and Ai, and as a result, the action functional within the perturbation theory 
in the deformation parameter A. The main result of this paper is the derivation of 
the modified Einstein equations in ( 16.41 ) in the weak deformation limit. In this case 
the pure noncommutative fields, namely /i'"' and cr, contribute to the right-hand side 
of the Einstein equation, that is, the stress-energy tensor The explicit form of these 
non-commutative correction terms has been derived in ( 16.71 ) through ( 16.121 ) for the first 
time. 

It is worth noting that the dynamical equations for the fields, found in this paper, 
are classical and therefore they should be studied from the classical point of view. It 
is an intriguing question to study the physical and mathematical effects of the non- 
commutative corrections. First of all, this is the problem of singularities in general rel- 
ativity. Another very interesting task would be to find some basic particular solutions, 
in particular, static spherically symmetric solutions (non-commutative Schwarzschild) 
and time-dependent homogeneous solutions (non-commutative Robertson- Walker). Fi- 
nally, it is the question whether the non-commutative fields could account for the dark 
matter and the dark energy. All these questions require further study. 

Some of the physical implications of Matrix Gravity have been extensively dis- 
cussed in ISl l9l [Toll . This theory exhibits non-geodesic motion which can be related to 
a violation of the equivalence principle, moreover, because of the new gauge symmetry, 
there are new physical conserved charges. At last, this theory represents a consistent 
model of interacting spin-2 particles on curved space which usually was a problem. An 
interesting question is the limit as ^ oo of our model, this might be related to matrix 
models and string theory. 

As it is outlined in the introduction Matrix Gravity can be considered as a Grav- 
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itational Chromodynamics describing the gravitational interaction of a new degree of 
freedom that we call gravitational color. Whether or not it is related to the color of 
QCD is an open question. Let's suppose for simplicity that it is the same, and that 
the gauge group of Matrix Gravity is nothing but S U{3). If one pushes this analogy 
with QCD to its logical limit then this would mean that the theory predicts that the 
gravitational interaction of quarks depends on their colors. Exactly as in QCD the 
strong interaction between quark of color / and a quark of color j is transmitted by 
gluon of type (//), the gravitational interaction between quark of color / and a quark 
of color i is transmitted by the graviton of type {if). In this case, all particles in the 
electro-weak sector, including photon, do not feel the gravitational color. In that sense 
it is 'dark'. Notice that in the non-relativistic limit the Newtonian potential will also 
become 'matrix- valued' . One can go even further. Since the usual (white) mass is de- 
termined by the sum of the color masses, one can assume that the color masses can be 
even negative. Then the gravitational interaction of such particles would include non 
only attractive forces but also repellent forces (antigravity?). This feature could then 
solve the mistery of singularities in general relativity. 

The consequences of our model in the ambit of cosmology are easily seen by in- 
specting equations ( I6.4l i and ( I6.5l l. The deformation of the energy-momentum tensor 
in ( 16.51 ) is written only in terms of the noncommutative part h^^^ . It would be interesting 
to study whether or not /i'"' could account for a field of negative pressure. If so, our 
model could describe the dynamics of dark energy. Furthermore, the distortion of the 
gravitation expansion of the universe due to the non-commutative degrees of freedom 
of the gravitational field will certainly have some effect on the anisotropy of the cosmic 
background radiation, nucleosynthesis and structure formation. Of course, a detailed 
analysis of these effects requires a careful study of the fluctuations in the early universe. 

Of course the validity of this statements require further investigations. The ultimate 
goal of this theory is to construct a consistent theory of the gravitational field which is 
compatible with the Standard Model and able to solve the current open issues which 
afflict General Relativity, i.e. the problems of the origin of dark matter and dark energy, 
the recent anomalies found in the solar system (Pioneer Anomaly, flyby anomaly, etc.) 
and last, but not least, the problem of quantization of the gravitational field. 

In summary, we would like to stress that our model makes it possible to make a 
number of very specific predictions that can serve as experimental tests of the theory. 
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Appendix A. Coincidence Limits 

In this appendix we will list the various coincidence limits used during the calculations 
performed in this paper. 

Through all this appendix the subscripts 0,1 and 2 will be used to denote terms 
of different order in the deformation parameter A. More precisely for any quantity X 
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which contains different orders of A we write 

where A'o, X\ and X2 denote, respectively, the zeroth, first and second order in A. 

We start with the coincidence limit of the operator H in ( 13.161 ) and its derivatives. 
More precisely we have 

m^^^.^vh"" . (A.l) 

For the first derivative we obtain 

[H;,]^A^,4ph''P,^. (A.2) 
For the second derivative we get the following formula 

[H;^,A = -^^,^PR%yp - jA^a^h'-^iR^yp + R%vt;) + A^a^ph^P.^y . (A.3) 

Recall, now, the definition ( 13.171 ) for the operator K. The coincidence limits of the 
terms in K are 

= -2gP\ (A.4) 

[b^'"] = -2hP\ (A.5) 

For the derivatives of these quantities we have 

[(V;,B''''')o] = 0, (A.6) 

[(V^B'''Oi] = -2//";^. (A.7) 

For the other terms we have 

[g'o] = 0, (A.8) 

[g{] = -hP'',,. (A.9) 

for the derivatives of (IA.8I 1 we get 

[(Vv.G^')0] = ^RPy+liPy, (A. 10) 

[(VyGP'h] = -hP''.^^-ihP''R„y+^hP-'RP,y,+hP''Ky , (A. 11) 

[(V.G^')2] = -[cr;,v,/0 - " [cr''^;„cr] - [cr-P,(r,y] . (A.12) 

For the operator X. in ( 13.181 1 we need the following coincidence limits 

my)0] = -[(^v)0] =-f V, (A. 13) 

[(^v)0] = -[(^v)0] =-(^;v, (A. 14) 

my)o] = my)o]^\[cr;y,cr]. (A. 15) 
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We also used, during the calculation, the coincidence limits for the derivatives of Jiy, 
namely 

[(V^^v)ol = -<P:vti+^^Rvti-\'Ryi,, (A. 16) 

[{^^,^y)l\ = -(Ty,, (A.17) 
[(V;.^v)2] = \{[(r;r„cT] + [cT.,y,a.,,]). (A.18) 

Appendix B. Integrals 

In this appendix we will explicitly derive the solution to the integrals in ( 14.35b and 
(14.421) . The first integral that we want to evaluate is the following 

1 T2 

/,(a,^,r,5)= J dT2 J dT,Tl(l-T2f{T2-Tiy(l-T2+T,f . (B.l) 


This integral is well defined for ReCa) > -1, Re(y) > -1 and < T2 < 1. By using the 
integral representation of the hypergeometric function we can evaluate the integral in 
Ti, i.e. 

1 

r(i + g)r(i + y) r i+^+y,, -js+s 1. ^ ^2 \ 

r(2 + a + y) J \ 1 - T2 / 



(B.2) 

Now, by using the linear transformation formula for the hypergeometric function we 
get HI 

/ T2 

2F1 \ \ + a, —5, 2 + a + y ; 



1-T2 



= (1 - T2)'+"2/^i (1 + a,2 + a + y + 6,2 + a + y ■T2) . (B.3) 



By substituting the last expression in the integral (IB.2l i we obtain 
r(i + a)r(l +y) 



r(2 + ff + y) 

1 



^ dT2T^"^^{l - T2)^^''^'^^'\Fi (1 + a,2 + a + y + S,2 + a + y ■,T2) . 



(B.4) 

This integral over T2 is, now, of the following general form 

1 

1 = J dx x''-\l - xy-^2Fi(a,b,c -x) , (B.5) 
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which is well define for Re(c) > and Re(t/) > and has a solution in a closed form 
ifTSl . namely 

nc)nd)T{c + d-a-b) 
T{c + d-a)Y(c + d-b) ' 
By using the results ( IB.6I ) in the integral ( IB.4I ). we get the solution ( |4.36t . i.e. 

J ^ r(l + a)Y{\ +/j)r(l + y)r(2 + a+li + 5) 
' r(3 + ff + j6 + y + S)Y{2 + a + jS) 

The next integral that we used in our paper is the one in (14.421 1. namely 

1 

/2(a,y6) = J dT,T",(\ - T,f . (B.8) 



The solution to this integral is easily find by recalling the integral representation of the 
hypergeometric function [[T] 

1 

2^1 {a,b,c;z) = ^^J^f r dt t'-Hl - tr-'-\l - tzT" , (B.9) 

r(fo)r(c -b) J 



where Re(c) > and Re(b) > 0. From this last general expression we obtain the 
integral (IB. 8b by setting z = Q, a-b-l and [i - c - b - I. By recalling that 
2Fy {a,b,c ;Q) - I, we finally get 

^' r(2 + «+/?) ^ ^ \ I 

where B(a, b) denotes the Euler beta function. 
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